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0\ ; Abstract 



We employ the Schwinger boson mean-field approach to study the effects of 
arbitrary frustrated bonds and plaquettes (formed from four frustrated bonds) 
in two-dimensional ferro- and antiferromagnets on the spin-wave spectrum 
and the correlation length at finite temperatures. We distinguish between 
strongly frustrated bonds (plaquettes), when the frustrated coupling J' ex- 
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ceeds the spin canting threshold J c , and weakly frustrated bonds (plaquettes) 



with J' < J c , (J c — J')/J c ~ 1. It is shown that in antiferromagnets the ampli- 

0' 

tude of spin-wave scattering on strongly frustrated bonds or plaquettes grows 



with the decrease of the temperature. A small amount of such defects re- 
duces significantly the spin-wave stiffness and the correlation length at low 



/\ ' temperatures. As a result, the quasi-2D Neel temperature is sharply sup- 



pressed. Quantum fluctuations are also considered and their effect on the 
spin- wave spectrum is shown to be of the order of (2S)~ 2 ln~ 2S in the large 
spin limit. For weakly frustrated (nonfrustrated) defect bonds (plaquettes) 
the spin-wave stiffness renormalization is of the order of the dopant concen- 
tration and does not depend on the temperature. The results account for the 
observed properties of doped quasi-2D La2Cu04 +x . The frustrated bonds in 
2D ferromagnets act in the same way as in antiferromagnets, while the effect 



of frustrated plaquettes is weak and temperature independent. 
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I. INTRODUCTION 

The discovery of high T c superconductivity in the doped cuprates La 2 CuO i+x and 
YBa 2 Cu^O^ +x enhanced the interest in the study of two-dimensional (2D) quantum an- 
tiferromagnets. It is now well established that the system of localized Cu 2+ spins in both 
families (near x— 0) is described by the 2D Heisenberg Hamiltonian, and a small coupling 
between Cu 2+ layers leads to the formation of a three-dimensional Neel ordered state at 
temperatures much smaller than the intraplane exchange interaction energy.a'B Above the 
Neel temperature T/v one observes 2D antiferromagnetic correlations. To leading order, the 
correlation length £ in undoped samples followsall the theoretically predictedoB law 

£cxexp(A/T). (1.1) 

The calculation by Chakravarty et al.p based on the quantum nonlinear sigma model 
(QNLaM), showed that A = 2irp s , where p s is the spin stiffness constant. The preex- 
ponentional factor in the two-loop approximation does not depend on the temperature up 
to terms of the order of T/2irp s in antiferromagnets and scales as T 1//2 in ferromagnets. A 
very good quantitative agreement between the QNLaM correlation length and the results 
of a neutron scattering study of the model 2D, S = 1/2, square-lattice Heisenberg anti- 
ferromagnet Sr 2 Cu0 2 Cl 2 was recently obtained in a surprisingly wide temperature range 
T N <T < 2tt Ps B 

Another approach to low- dimensional (D = 1,2) magnets has been proposed by Arovas 
and AuerbachB They employed the Schwinger boson representation for spins in a SU(N) 
generalization of the Heisenberg model. Their expression for the correlation length corre- 
sponds to the one-loop aproximation in the QNLaM. The constant A in Eq. (|1.1|) is 27rJS 2 
in the limit of large S, and J stands for the intraplane coupling constant. 

The doping has a drastic effect on the magnetic properties of the cuprates. Even a very 
small dopant concentration, which indroduces holes residing on the oxygen orbitals in the 
Cu0 2 planes, substantially reduces the Neel temperature JiBBE! In La 2 Cu0 4 , doped with 



strontium or with excess oxygen, the long range order disappears at doping concentration as 
small as 2%. It has been shownH'Elu3 that in samples of La,2CuO± +x with a small concentration 
of defects (of the order or less than 1%) the correlation length varies little with doping at very 
high temperatures 500 - 600 K, however the growth of the spin correlations is increasingly 
inhibited by the defects as the temperature is decreased. 

In this paper we shall focus on the properties of doped 2D magnets. Doping with excess 
oxygen or with strontium creates holes on the oxygen sites in the Cu02 planes.t!-!J Since the 
Cu-0 distance is half the Cu-Cu distance, the exchange interaction, J CT , of the hole spin with 
the two nearest-neighbor Cu spins is much larger than the Cu-Cu exchange, J. For either 
sign of J a the hole spin would thus like to be parallel (antiparallel) to the two neighboring Cu 
spins. Therefore the oxygen hole introduces an effective ferromagnetic coupling, J', between 
its two Cu nighbors, with J' = 0(\ J a |) 3> J ■ This was the basic idea of the frustration 
model of Aharony et al.x-i 

Motivated by these arguments and the results of the experimental investigation of 
La2CuO^ +x and La 2 - x Sr X C 'uO ' 4 pntl'&3'El we shall consider two simple models of defects, 
namely, frustrated bonds and frustrated plaquettes. We hope that the main qualitative 
results, presented in this paper, survive in the case of more complicated frustrating defects. 

The effect of randomly distributed frustrated bonds on the Heisenberg ferro- and anti- 
ferromagnets in the ordered state has been investigated by many author s.EJlia It has been 
shown that there is a threshold energy of the frustrated bond, 

J c = (z/2-l)J, (1.2) 

which determines the local stability of the system. Here z is the number of the nearest 
neighbors. For a 2D square lattice, z = 4 and J c = J. If the (positive) energy of the 
frustrated bond J' exceeds J c , the two spins connected by the frustrated bond cant. In the 
classical limit the canted spins act on the magnetic background as a dipole.li^Hi^l As a result, 
the spins around the ferromagnetic bond are also canted, and the canting angle in a 2D 
Heisenberg magnet decays at large distances r from the defect as r _1 . 



The quantum fluctuations of the spins, connected by the frustrated bond, have been 
calculated for J' < J c in the framework of the linear spin wave theoryJiZnla These fluctuations 
diverge when J' —>■ J~ , which reflects the breakdown of the linear spin wave approach. 

The effect of localized holes on the properties of quasi-2D antiferromagnets at finite 
temperatures has been studied by Glazman and Ioselevich.By A classical approach, based 
on the dipole model of Aharony et al.f^ led them to the conclusion that the renormalized 
stiffness is a function of (x/T), where x is the concentration of localized holes. 

The model of a single frustrated bond is reasonable for samples doped with excess oxygen. 
In contrast, Sr doping replaces a trivalent La 3+ ion with a divalent Sr 2+ ion in the plane 
above a C11O2 plane. The Coulomb potential that pins the hole in the CuOi plane is 
then due to the Sr 2+ ions, which project onto the centers of the copper-oxygen plaquettes. 
Therefore the holes are localized on small regions around the centers of the above plaquettes. 
The simplest defect model of this kind is a plaquette formed by four frustrated bonds£lrc3 
The canted spins at the corners of the plaquette act on the long-range order parameter (at 
T — > 0) as a quadrupole rather than a dipole. The energy spectrum of the hole, localized on 
a plaquette, was calculated in Refs. |2TJ and |22|. To our best knowledge neither the critical 
value J c of the frustrated coupling (at T=0), nor the effect of the frustrated plaquette at 
finite temperatures, has been considered till now. 

In this paper we study the effects of defect (frustrated and nonfrustrated) bonds and 
plaquettes in 2D ferro- and antiferromagnets on the spin wave spectrum and on the corre- 
lation length at finite temperatures. The Schwinger boson mean-field approach (SBMFA)Q 
will be used. This method has been applied successfully to study a variety of properties of 
low dimensional antiferromagnets.c3'B3i^l 

We show that there is a striking distinction between the effect of strongly and weakly 
frustrated defects on the behavior of doped antiferromagnets. For strongly frustrated defects, 
with t = (J'—J c )/Jc > (J c = J for frustrated bonds, J c = 0.376 J for frustrated plaquettes), 
the renormalization of the spin wave velocity c and the correlation length in the large spin 
limit (to leading linear order in x) is given by 



c(x) - c(0) AJSH 

~^T~ = - a WT7f> (L3) 

ln^Il = -a 87lj2SH X (14) 

111 £(*,0) V(l + t) X ' li - 4j 

where the multiplier a depends on the type of the defect (a ~ 1 for frustrated bonds, a ~ 
3.5 for frustrated plaquettes), and decreases slightly with the decrease of the temperature. 
Hence, a small amount of frustrated bonds or plaquettes reduces significantly the spin-wave 
velocity and the correlation length at low temperatures T <C AJS 2 , if t is of the order or 
larger than unity. The decrease of the 2D correlation length with doping causes the strong 
supression of the quasi-2D long range magnetic order. The Neel temperature changes with 
the dopant concentration as 

T N (x)-T N (0) US 2 

— ' , . — — ~ —a — , r x. (1.5) 

T N (0) T N (0) K ' 

Because of the large factor 4JS 2 a/T N (0), the Neel temperature extrapolates to zero at a 
doping concentration which is much smaller than the percolation threshold. 

The quantum fluctuations, associated with the frustrated bonds, are relatively large 
near the threshold J' = J c . However, when J' ^> J c , they are surprisingly small, of the 
order of (2,S)~ 2 ln~ 2S. For weakly frustrated (nonfrustrated) defect bonds or plaquettes, 
with (J c — J')/J c ~ 1, the spin-wave stiffness renormalization is of the order of the dopant 
concentration and does not depend on the temperature. 

The frustrated bonds in 2D ferromagnets act in the same way as in antiferromagnets. 
However the effect of frustrated plaquettes in ferromagnets is entirely different. The renor- 
malization of the spin-wave stiffness and, hence, of the correlation length is not enhanced, 
and it is temperature independent at both J' < J c and J' > J c . 

The above results account for the observed properties of doped quasi-2D antiferromag- 
nets. 

The paper is organized as follows. In Sec. |TI| the SBMFA is used to study the problem of 
a defect bond in a 2D Heisenberg ferromagnet. We derive the mean-field Hamiltonian and 



find the temperature dependence of the spin-wave scattering amplitude on the defect bond. 
Then we treat the ferromagnet with a low concentration of bond defects, and obtain the 
renormalization of the spin-wave stiffness and the correlation length at small temperatures 
T <C 2irp s . In Sec. |J.| we extend this method to treat 2D antiferromagnets with defect 
bonds, and study the effect of such defects on T N . The plaquette type defects are considered 
in Sec. |IV|. In Sec. [V| we summarize our results, and compare them with the experimental 
findings. Some comments are made concerning the temperature and defect concentration 
dependence of the spin wave velocities in quasi-2D antiferromagnets at T <C T/y. The canting 
of the spins at large distances from the strongly frustrated (J' > J c ) bond is addressed in 
Appendix A. In Appendix B we calculate the T-matrix for the scattering of the Holstein- 
Primakoff spin waves on frustrated plaquettes at T = and J' < J c . 

II. FRUSTRATED BONDS IN A 2D FERROMAGNET 

A. The Hamiltonian 

We consider a Heisenberg ferromagnet on a square lattice with nearest neighbor (nn) 
interactions only. Thus, in the standard Heisenberg Hamiltonian 

H = — 2_j JijSi ■ Sj, (2.1) 

the sum is over nn pairs, < ij >, and the exchange interactions J^- are equal to J (for most 
bonds) or —J' (for a small concentration of bonds). 

In the Schwinger representation, each spin is replaced by two bosons: 

S+ = a\b u Sr = ai bl S? = I(aIa i -&J& i ), (2.2) 



together with the constraints 



^{a\a i + b\b i )=S. (2.3) 



After this transformation is performed, the Hamiltonian can be written as3 



H=-2J2 Jij ■ BtjBij : + J2 Hahi + b\h), (2.4) 

where B\, = \{a\aj + b[bj), and : : denotes normal ordering. The second term in this 
expression is added to take into account the constraint (|2.3|) , Aj are the Lagrange multipliers. 



An important approximation in the mean field theory of Arovas and AuerbachQ is that the 
constraints Q2.3|) are imposed on the average, i.e. they are taken into account by introducing 
a single Lagrange multiplier A. This approximation is plausible in a spatially ordered magnet, 
but in a disordered one the dependence of Aj on the position may be important. In what 
follows we introduce two Lagrange multipliers: Ai for the sites connected by the defect bond, 
and A for all other sites. 

Performing the mean- field decouplings in Eq. fl2.4|), one obtains 



H = H + H mt . (2.5) 

Here H Q is the Hamiltonian of the undoped ferromagnet, 

H = -2Q f J £ (B tj + 4) + A ]T(aJa, + b\b t ), (2.6) 






where the mean-field amplitude Qf =< B], >=< B^ > describes the short-ranged ferromag- 
netic correlations.tBEa The amplitude Qf, like A, is supposed to be positionally independet 
in the ordered crystal B 

The term H int in Eq. ( |2.5| ) gives the interaction of the spin waves with the defect bonds, 
< Im >, 

H int = 2W f £ (B lm + Bjj + (A x - A) £ (a] ai + ala m + b\b L + bl n b m ), (2.7) 

<lm> <lm> 

where 

W f = Q f J + Q 1 J'. (2.8) 

The mean-field amplitude Q\ =< B^^ l+ ^ > describes the ferromagnetic correlations of 
the spins at the ends of the frustrated bond a. It differs from the amplitude Qf for the 
unfrustrated bonds, and this difference is crucial for the properties of the doped magnets. 
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Transforming Eq. ( |2.6|) to the momentum space, we find 

# = $>*( a W +6 W)- ( 2 - 9 ) 

g 

The excitation spectrum u^ has the form 

W3.= A + JQ/z(l-7 f ), (2.10) 

where 7g> = z ^2^exp(iq • 5) = (cosg x + cosq y )/2, z = 4 is the number of nearest neigh- 
bors, and 5 is summed over nn vectors. The spectrum is characterized by two temperature 
dependent parameters, namely, the mean-field amplitude Qf and the gap A = A — JzQf. 
They are governed by the equations^ 

^E«K) = s (2.ii) 

Q/ = S-^£»(wf)(l-7*), (2-12) 



which follow from the constraint (|2.3| ) and the above definition of Qf. Here n(u$) = 
[exp(u)f/T) — 1] _1 is the Bose distribution function, N is the total number of spins. 

It is seen from the last equation that Qf is equal to S up to corrections of order (T/AJS) 2 . 
Then Eq. (|2.11| ) yields that the gap A is finite and exponentially small at low temperatures, 

A oc T exp(-4vr JS 2 /T). (2.13) 



The formula (|1 . 1|) for the correlation length £(T) follows immediately from this expression 
and Eq. (pT0|).l 

The 2x2 interaction matrix for each defect bond in Eq. (2.7) can be easily diagonalized. 
The interaction Hamiltonian (2.7) is then rewritten as 

H mt = {2JQ f Hf - W f ) £ [(oj - al)( ai - a m ) + {b\ - 6+ a )(6, - b m )\ 

<lm> 

+2JQ f fi f J2 [(«! + fl L) (ai + Am) + (&J + &L) (&l + &m)] , (2. 14) 

</m> 



where /i/ = (Ai — \+Wf)/4JQf. Performing the Fourier-transformation of this Hamiltonian 
we find 

Him = E v v E ^ E e <( *-* ) ^/,(fi)/;(fi)(a^a & + &}.&&)■ (2.15) 

Here 

A(g) = 1 - e"^, / 2 (g) = 1 + e~^ s , (2.16) 

i>i = 2JQf/j,f — Wf, v 2 = 2JQf/j,f, and I is summed over the defect bonds. 

B. The single defect problem 

In this subsection we use the Green's function method to solve the problem of spin wave 
scattering on a single defect bond. Let us define the retarded Green's function by the relation 

G(q,&,t-h) = -ie(t-ti) < Mt),<4(*i)l > T (2.17) 

Here ©(£) is the step-function and < • • • >t denotes the thermodynamic average. 

In a ferromagnet with one defect bond, the Fourier-transform G(q, ft, uj) is linked to the 
T-matrix by the relation 

G(ftft,u;) = G°(q,w)8 m +G (q,u)T(q,q 1 ,u)G°{q u w), (2.18) 

where G°(q,u) is the Green's function for the pure ferromagnet, 

G (q,Lu) = (iu-LU (! )-\ (2.19) 

and the T-matrix satisfies the equation 

T(ft ft, u;) = V(q, ft) + 1 E V(ft g 2 )G°(g 2 , u;)T(ft, ft, u), (2.20) 

92 

with the interaction energy 

V(q, ft) = E vJMfM)- (2.21) 



i^=i 
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The Green's function for the 6-operators obeys the same equations (|2.18 ) and Q2.20| ). 



The constraint ( |2.3| ) for the spins, at the ends of the frustrated bond, can be written in 
the form 



\ /-+oo \ 

<a]ai> T = / n(u>)— J^ ImG(q 1 ,q 2 ,uj)duJ = S. (2.22) 

77 J-oo iV t~Z 



l r +oa , , l 

9192 

We put here for sake of simplicity r\ = 0. Substituting into this equation the Green's 
function from Eq. Q2.18|) , and taking into account that 



/+oo 1 1 

n^Im- £ G°(q, u) = - £ n(w f ) = S, (2.23) 



we find 



/+oo 
dun(ij)Y / ImG (q,ij)T(q,q 1 ,uj)G (q 1 ,uj) =0. (2.24) 

-00 =" 



r+oo 

dun(u) 
In the same way, from the definition of Q\ we get the second equation 



Q 1 = < a \a m >= Q f 

1 /-foo 1 

+ - / dwn(u)— Y l ImG°(q,u;)T(q,q 1 ;u;)G (q 1 ,u;)(l-e^ a )- (2-25) 

991 

When deriving this equation the constraint ( |2.24|) has been taken into account. 

We now solve the set of Eqs. (|2.20|) , ( |2.24j) and (|2.25|) . The perturbation energy V(q, q\) 



is the sum of separable terms, and therefore the integral equation (|2.20|) is solved easily. 
The result is 

Tm l )=Y. v " f f!'t ) - < 2 - 26 > 

Here 

D v (u) = l-v v <f> v {u), (2.27) 

where 

M«>) = ^Y.fMG°(q,uj)f „(<!)■ (2.28) 

q 
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The functions 4> v {oj) can be rewritten with the aid of Eqs. Q2.16|) , ( |2.19| ) in the form 



0i (w) = ^j^[-l + U(p{u)], 



MUJ) = 2JS + {A -2JS MUJ) - 



(2.29) 



Here 



v(") = ^ 5> - "aT 1 (2-30) 



is the local Green's function for the undoped ferromagnet. We find from ( |2.30| ) that at small 
uj (to < AJS) 

%M = '4^' V 4 -A| ' /Tr^(w) = -^e(o;-A). (2.31) 

Eqs. (|2T2l) - (P^0|) ) yield 

^(w) = C/(T) + /i/ + [1 - C/(T) - M/W(u>), (2.32) 

where 

WV(T) 

C '< T) = ' - w- (2 ' 33) 

We show next that ///(T) is always small at small temperatures T <C 4JS 1 , and the function 
Cf(T) is small, if the frustration is sufficiently strong. Thus, the function D\{uS) is small 
too, when the defect bond is frustrated. This means, that the T-matrix and, hence, the 
spin wave scattering amplitude are enhanced strongly in this case. This is the reason of the 
unusually large effect of frustrated bonds on the properties of the doped magnets. 

We substitute now the function T(q\, q?) from ( |2.26|) into Eqs. ( |2.24| ) and ( |2.25| ). Taking 
into account that 

±J2G°(q,u)f n (q) = ^(w), (2.34) 

q 

we rewrite these equations in terms of the functions 4> u (lo) 
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t>i / dcun(uj)Im 



DAuj) J-oo K J D 2 (uj) 



Qi = S + — W dun(uj)Im—-^- 

271 J-oo DAi. 



2 
L 



^2.35) 



^2.36) 



When writing Eq. (|2.36| ) we have replaced the parameter Qf for the pure crystal by S, since, 
as we shall see later, the difference S — Qf ~ (T/AJS) 2 is always smaller, than the second 
term in the r.h.s. of this equation. 

It follows from Eq. (|2733|) and the definition (gj) of W f that 

-1-2C/ 



Qi = S- 



t + l 



(2.37) 



where t=(J'-J)/J. Eqs. (F36[) and (|2~3l) yield 

t + 2C/ " " +no --- 2 



t + 1 






(2.38) 



After the expressions ( 2.29|) - Q2.31 ) for the functions 0i(o;) and 02 (<^) are inserted into Eqs. 
Q ) and (^38|) , one finds 

2JS 2 /i/ 1 



T + 87^/1/ 85 



tf(T) 



* + 2Q 1 K(T) 



(2.39a) 



(2.39b) 



where 



x , 1 



2 i ™2l-l 



("OO ry, 

K{T) = - / dxxn(USx) [{( + fit - - In -) 2 + rr 

7T JO ' 7T X 



^2.40) 



When deriving these equestions we neglected terms of higher order in the small quantities 
Q and /i/. 

In the quasiclassical limit, S 3> 1, we find that with logarithmic accuracy 

Tl(T) 



K(T) 



8JS\Cf(T)+ti f {T) 



(2.41) 



where 



13 



2 11 
l(T) = 1 + - arctanf- In ], 1 < l(T) < 2, 

7T 7T C/ + tV 



if (f + fif > 0, and 



(2.42) 



2 1 1 

l(T) = 1 arctanf— In — r l, 



if C/ + /// < 0. 

From Eqs. ( [2.39|) - (|2.41|) we finally obtain 



. l(T), T , 2/1 8J5%. 



^2.43) 



(2.44a) 






(2.44b) 



Eqs. (|2.44j ) have a solution with (f + /i/ > at both £ > and t < 0: 



C ' (T) = J<"* + wsfei 



+p+ ™<2£±v } , 



(2.45a) 



/'/ 



8JS 2 t + 2 ' 



(2.45b) 



and 



A(0) = C/ + H = \{~t + [t 2 + ^(t + l)] 1 / 2 }- 



^2.46) 



If J' < J (t < 0) Eqs. ( F57D and (|2T45|) yield 
1)T/JS 2 > t 2 







Tl ( T ) \l/2 + Qiyv^Lw 2 ! 



16 JS 2 ' 



T 






(2.47a) 



ft - *[! " (^) 1/2 1 



(2.47b) 



2) T/JS 2 < t 2 < 1 







t I TZ(T) 



+ 0( 



r- 



2 8JS 2 It I V J 2 5 4 It I 3 



/'.f 



T 2 /(T) 



(8JS 2 ) 2 | t 



«C/ 



(2.48a) 
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Ql = s{1 - 4J§^T\ ] - (2 ' 48b) 

The function Q\ decreases with the increase of T much faster than Q. When T — > 0, 
-Di(O) = (f =\ t \ /2, and Q±, like Q, tends to S. Thus, the correlations of the spins, 
connected by the frustrated bond, are purely ferromagnetic, when the frustrated exchange 
J' is lower than the threshold J c = J. 

The properties of the defect change drastically, when J' > J (t > 0). If T j .IS 2 <C t 2 , we 
obtain 



It follows from Eqs. <^M) and ( ^49|) that Qi(T -»■ 0) = S/(t + 1) < S. Thus, when 



J' > J, the ferromagnetic correlations decrease with the increase of J'. At large J' 3> J, 
the correlation amplitude Qi(0) -C S, i.e. the spins are almost antiparallel. The correlation 
function Q^v +( f of two neighboring spins at large distances r from the frustrated bond is also 
less than S, and approaches S with the increase of r as r -4 (See Appendix A). Thus, the 
spins at the ends of the frustrated bond act on the ferromagnetic background like a dipole, 
in agreement with the classical picture.^ 4 ! 
The negative solution of Eqs. (|2.44|) 



Q = D x (0) = -t/2 + 0(T/US 2 ), t < 1, (2.50) 



which appears, when t > 0, resembles the bound state solution, obtained in Ref. 1 4 for a 



frustrated bond in a 3D ferromagnet. Indeed, the T-matrix ( |2.26| ) in this case diverges at the 
negative energy uj = —2irJSt/ln(l/t), which is the condition of a bound state formation.E3 
At t ~ 1 this state lies far from the bottom of the spin-wave band. Hence, it almost does 
not alter the spin- wave spectrum. In contrast, the dipole-type state affects strongly the 
properties of the ferromagnetliHHi3. Therefore we consider in what follows only the dipole- 
type states. 
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Eqs. (|2.39|) can be solved also at arbitrary spins. If | t |<C 1, the function /ij <C Q, and 
(f = -Di(O) is governed by the equation: 



C / 2 (2C / + t) = 0.5245(J^) 2 . (2.51) 



This equation yields 

1)J' <J, | t |> (T/AJS) 2 / 3 



« T ' = T + 1 4 ( ^ + 0^ (2 ' 52) 



2) | t |< (T/4JS) 2 / 3 



C/(^) = ^(7^) 2/3 + 0(M), (2.53) 



3) J' > J, 1 > t > (T/4JS) 2 / 3 

0.7245 1 / 2 T „ r , T 



« T » = ^i^ + °iW ! i' (2 - 54) 



Finally, if £ > 1, the function D± is given by Eq. ( |2.49b|) with l(T) = 2, up to small terms of 



the order of ln~ AJS/T. Thus, at t > the main features of the functions C/CO an d Qi{T) 
are the same, as for large spins. In other words, the dipole-like picture of the frustrated 
bond defect is valid not only in the classical limit, but rather at arbitrary spins. 

C. The renormalization of the spin- wave spectrum 

Averaging Eq. (|2.18|) over the distribution of the defects and the orientations of the 
defect bonds, we find the renormalized spin wave spectrum e$. To first order in the defect 
concentration x the configurationally averaged Green's function is given as 

G _1 (g, uj) = uj — ujq— S(g, to), (2.55) 

where the self-energy is given by E(g, u) = xT(q,q;uj). Thus, 

eg-= u^+xT(q,q;u^). (2.56) 
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In the long- wave limit (e^ <C T/ In (4 JS/e^)), we obtain for e^ 

e?=n q il-^-] + A 1 (T), (2.57) 

where Qq- = 4 J5(l — 7^). The properties of the renormalized gap A x will be discussed later. 
When the frustration is weak, we obtain -Di(O) ~| t |« 1. The renormalization of the 
spin- wave stiffness is of the order of x and temperature independent. At strong frustrations 
one has -Di(O) ~ C/CO ^ 1- The renormalization is enhanced, and the enhancement 
increases with the decrease of the temperature (see Eqs. (|2.47| ) - (|2.49|) , (|2.52|) - ( [2.54|) ). 



It follows from Eqs. ( [2.57|) and (|2.49b|) that for strongly frustrated bonds (J 7 > J) the 



renormalization of the stiffness p s can be written as 

P^-P^ = %, (2.58) 

where Uf = 8JS 2 t/l(T)(t + 1). Even a small defect concentration x < T/Uf <C 1 reduces 
significantly the stiffness. 

D. The correlation length 

The gap in the spin wave spectrum of a pure ferromagnet is governed, as mentioned above, 
by the constraint (|2~3l ). A self-consistent way to obtain the gap in a doped configurationally 
averaged crystal is to impose on the spectrum ( p. 57 ) the constraint (|2.3| ), averaged over the 



defect distribution. Since the defects under consideration preserve the value of the spins, we 
obtain the same Eqs. (|2.11|) , (|2.12|) as for the undoped crystal, with u^ replaced by e^. The 



gap Ax is then given by Eq fl2.13|) with the exchange energy J replaced by the renormalized 
value J(l — x/Di(0)). Thus, the correlation length in the doped crystal is given by 

£(T,*)cxexp[^^(l-— )]. (2.59) 

At a given x the ratio £(T, x)/C,(T, 0) decreases rapidly with the decrease of the temperature. 
Even though the expression fl2.57| ) is valid only at small concentrations x <C -Di(O), the renor- 
malization of the correlation length may be exponentially large, if x ^> D 1 (0)(T/27rJS 2 ). 
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III. DEFECT BONDS IN ANTIFERROMAGNETS 

A. The Hamiltonian 

We shall consider a bipartite antiferromagnet on a square lattice. Only nearest neighbor 
interactions are included in the Hamiltonian, and it is supposed, as before, that there is a 
small number of defect bonds with exchange energy J' ^ J, where J is the host exchange 
energy. The Hamiltonian is given by 

H=Y,JiA-Si, (3-1) 

where the couplings Jy are equal to J or —J'. The transformation to Schwinger bosons for 
the spins on the sublattices A and B is as follows: 

S\ = a^b, Sj = ab\ S Z A = -{a^a - tfb) (3.2a) 

S+ = -ab\ S B = -a\ S Z B = -{b ] b - a ] a) (3.2b) 

with the constraint ( [2.3|) . After the mean-field decouplings is performed, the pure antifer- 
romagnet Hamiltonian H and the interaction Hamiltonian H int transform to 

H = — — J2 ( a i a j + a \ a ) + b i b j + h \ h )) 

+ \J2(4^ + b\b t ), (3.3) 



H mt = (W- 2JQfi) ^ [(o, - al)(a m - of) + (6, - 6j n )(fo m - tfJ 

<lm> 

+2JQfi Y, [(«f + « m )( G m + ai) + (6} + b m )(b ] m + k)} . (3.4) 

<lm> 

Here W = JQ + J'Q', \i = (1/AJQ)(X 1 -\ + W), 

Q = - < didj + b^j >= - < ajo] + blbj >, (3.5) 



and Q' is given by Eq. ( |3.5| ) with the bonds < ij > replaced by the defect bond < Im >. 
The sums in (B.4I) are over all defect bonds < Im >. 



In momentum space the Hamiltonian Hq can be diagonalized by the standard Bogoliubov 
transformation 

a^ = u^ag- + v ^al- , (3.6a) 

b$ = u$/3$ + v$]_q , (3.6b) 

where 

2 JQz + A + Wf 2 JQZ + X-UJJ 

u* = -, v* — -. 3.7 

q 2uf ' q 200$ V ; 

The quasiparticle energy has the form 

^ 2 = (JzQ) 2 (l- 7a 2 ) + A 2 , (3.8) 



where the gap is given by A 2 = (A + JQz) 2 — (JQz) 2 . From equations (|273| ) and (|3T5D-([3~8"|) 
the self-consistent equations, obeyed by the mean-field amplitude Q and the gap A, can be 
obtainedEl 



q 






(3.9) 



In the limit of large spin S ^> 1, one has the amplitude Q = S(l + 0.158/2S), the gap 

A otTexp(-27vJQS*/T), (3.10) 

and the correlation length 

f (T) oc exp(27r JQS*/T), (3.11) 

where S* = 5(1-0. 197/5). @ 

The interaction Hamiltonian Q3.4| ) can be written in the momentum representation in a 
simple form, 
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1 l 

H in t = -T7 E V » E < A g I Vu{q) >< VM) I Aft > . 



(3.12) 



Here we introduced the two-component column vectors | A$ >, \ r] v {q) >, and row vectors 
<A$\, < T)„(q) | as: 

( „ \ 



Aj>-- 



, <^-i=(4 



\ a -<?J 



q- , a -q) j 



(3.13) 



Vu(q) > 



i ■ 



and V x = 2 JQ/i - W, V 2 = 2 JQ/x. 



, <^(g)|=(l, (-l)"e 



= ig-a 



(3.14) 



B. The one-bond problem 



Let us define the Green's function matrix by the relation 



G af }{q, qi;t-h) = -ie{t-t 1 )< [A a? {t), AL (*i)] > T 



where A]* = a*, A 2 ff = a)-,. The unperturbed Green's function matrix is as follows: 



(3.15) 



G°(g» 



4JQ + LU AJQ 1<r * 



u 2 -ul 



\ 



(3.16) 



/ 



4 JQy$ AJQ - uj 

The one-defect problem is solved in the same way, as for the ferromagnet. The Green's 
function matrix and the T- matrix obey Eqs. ( |2.18| ) and ( f2.20| ), with the functions G and T 
replaced by the 2x2 matrices G and T. The solution for T(q, gi; uS) is 



T(q,qi,u) = 2^ 



D v {u) 



(3.17) 



Here 



A,(w) = l - KVvM, 



(3.18) 



and 



20 



VV = ^E < ^® I £°(?» I »U$) > • (3-19) 

9 



Using Eqs. ( gig) and (|3l6|) we find 

1 



r,( ' ;;) ""2iQ [1 " (CL;2 " A2 ^ M] ' (3 - 2 ° a) 



^ M = 27g + (16jg " ^T^' (3 ' 20b) 



where the local unperturbed Green's function g{uS) is 



The function g{u) is easily calculated at small frequencies uj <C AJQ: 

1 8 JO 1 lu 2 

Re9(uj) = -8^JW [ln |^-A*|V2 " ^ + °W )l (3 ' 22a) 

J"#(w) = "T^v^ " A ) " ©(-w " A)]. (3.22b) 

From Eqs. ( p.20| ) and (|3.21|) we find the function Di(o;), 

^(w) = CCO + ix + (1 - C - A*)<A(w), (3.23) 

where 

C(T) = l-^. (3.24, 

We see that the spin-wave scattering on frustrated bonds is enhanced strongly, since, as we 
show below, ((T) and n(T) are small. 

Let us proceed to the derivation of the equations, which govern the function Q' and the 
local Lagrange multiplier \i. In terms of the Green's function matrix, the local constraint in 
the momentum representation can be written as 

dun(u)Im22 < Xt I G{q t qi,u) - G (g>)% \ \] >= 0, (3.25) 

-00 



<i<n 
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where | Xf > is the spin 1/2 spinor with spin up. 
Equations (|2TT8| ), (|3T^ ) and ( F^25|) yield 



dum(uj)Im ^ 14- 



L T „(u;)L I/T (^) 



v=i 



A,M 



0. 



(3.26) 



Here 



L T «>M = L„ T (u;) = —Yl < XT I G°(q,u) | r^(g) > . 



(3.27) 



The equation for Q' can be derived in the same way. The first step is to transform the initial 
expression 



QQl 



a W-gi : s 



iq-a 



(3.28) 



with the aid of Eq. (|3.17| ), and the constraint (|3.25| ). We find 

\ r+oo 

Q' = Q / dujn{uj) 

7T J-oo 






XT I G°(q,u) | r} u {q) >, 



(3.29) 



where < xi I is the spinor with spin down. Taking into account the relation 
< Xi e% ^ a ~ Xi 1= ~~ < Vi(q) I) we fi n d 



Q =Q + - (kun(uj)Im - , \ Im 

TV J-oo J^l{U>) 



;w). 



(3.30) 



The functions L u ^(u) can be written in terms of the functions ip v {uj) and g(w), as 

t t \ ^i( w ) , /- ^ 



(3.31a) 



L 2T H = (w + 8JQ)g(u) - -M")- 



(3.31b) 



Using the relation between Q' and £ 



Q' = Q 



1-2C 
1 + t 



(3.32) 



and Eqs. (|37TCf ) , O), ( ^26]) , and ftTgCD we finally obtain 
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^Te^ + T^ < 3 - 33a » 



' ^ 7T J ( T )' ( 3 ' 33b ) 



where 



t + 1 Q 



I(T) = - [°° dx[x 2 coth ^F - xC] [(C + /i - - In ^) 2 + x 4 ]- 1 - (3-34) 

n Jo 1 ir x z 

In the large spin limit the integral I(T) is given by 

^ = simh + ^(C + ,)KC + ,)- |1/2 + ° (1) + ° (,n " ?'■ C + " > ° (3 ' 35) 

The second term in the r.h.s. of this equation takes into account the effect of quantum 
fluctuations. The slowly varying function l(T) is given by Eq. (|2.42|) with C/> A 4 / replaced by 

Ca 4 - 

Several results follow immediately. First, we see that /i, like /if, tends to 0, when T — > 0. 
Near the threshold, when | t |<C 1, the quantum correction to the mean-field amplitude 
Q' exceeds the correction (of the order of l/S) to the function Q in the undoped samples. 
Indeed, if S~ 2 ^ 3 <C| A |<C 1, and T = 0, the mean-field amplitude becomes 

Q ' = Q -( \t\Mi/\t\/ ' 2 - < 336 > 

The quantum correction increases with the decrease of | t | as (| t | In | t |) _1//2 , in agreement 
with the result obtained in the linear spin wave theory.E-3 However, while in the linear spin- 
wave theory the quantum corrections diverge when 1 1 I— > 0,li3lla we find a finite value in this 



case 



Q ' =w - WW 1/3] - (337) 

It follows from Eqs. ( |3.33b| ) and (|3.35|) that the quantum corrections, associated with the 
spin-wave scattering on strongly frustrated bonds (t ^> S^ 2 ^ 3 ), are of the order of S~ 2 In" S, 
i.e. smaller than the usual l/S corrections: 
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c - { ~r ) ' 1 'wwor' (338) 



«" = TT7' 1 " <T> 2 (2wW < M9 > 



When deriving these equations we neglected terms of order unity in Eq. ( |3.35|) . The effect 



of these terms is to replace 2Q in Eqs. ( |3.38|) , ( p.39|) by 2Q + k, with k«1. Thus £ "C 1 for 



any value of the spin. This implies that the quantum effects are not very important even 
for real spins S = 1/2. 

Neglecting the quantum fluctuations, we find that at T = 0, t > one has C(T) = 0, 
and Q' = 5/(1 + i) < S. Like in the ferromagnet, the function S — Q^^.g decays as r -4 at 
large distances r from the frustrated bond (see Appendix A). Hence, the spins at the end of 
the frustrated bond act as a dipole, when J' > J. 

At finite temperatures Eqs. ( ft.33|) and (|3.35| ) yield in the large spin limit that at t ^> £ 



one has 

The function -Di(O) coincides, up to terms of the order of 1/S, with Eq. (|2.49|) for -Di(O) in 
a ferromagnet. The solution of Eqs. ( |3.33j )- p.35| ) in the quasiclassical limit at small t also 
coincides with the corresponding solutions [Eqs. ( |2.47[ ), ( |2.48[ )1 in a ferromagnet, and will 
not be given here. 

So far, we discussed only the positive solution D\ > of Eqs. (|3.33|) . As for the fer- 



romagnet, these equations also have a negative solution at J' > J. However, unlike the 
ferromagnet, this solution does not correspond to a bound state, since the zero temperature 
T-matrix in the antiferromagnet has no pole at negative u. 

C. The spin-wave spectrum and the correlation length 

To leading order in x, the configurationally averaged Green's function is given by 
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[1 - G (q;u)f(q,q;u)x]G(q,u) = G (q;uj). (3.41) 

The renormalized spin-wave spectrum, which follows from the solution of this equation, can 
be written in terms of the T(q, g)-matrix elements as 

4 = loI + x[AJQ(T n + T 22 ) + ^(T n - T 22 ) + 4JQ 7g -(T 12 + T 21 )]. (3.42) 

We have at small frequencies, u$ <C 4JQ, 

In the small q limit, when -Di(O) ~^> (co>|/87r J 2 Q 2 ) \n(4JQ/u)^), the spin-wave spectrum is 

4=4[(l--^] + A?CT), (3.44) 

where £7| = (4JQ) 2 (1 — 7|). At higher frequencies the function -Di(O) can be dropped in 
Eq. ( |3.43|) , and we find 



4 = 4 + Stt J 2 Q 2 a; In" 1 ^. (3.45) 

The wave-vector dependence of the last term in the r.h.s. of this equation is subtle. The 
spin-wave spectrum, hence, acquires a concentration dependent gap, the spin-wave velocity 
being the same as in the undoped antiferromagnet. 

Like for ferromagnets, the renolmalization in the small q region is of the order of x for 
weakly frustrated bonds, and is enhanced and temperature dependent, when t < 0, | t |<C 1 
or t > 0. For strongly frustrated bonds (t > 0) we have 

e}=El(l-^)+Al(T), (3.46) 

where 

4JQ 2 £ 

The renormalized spin wave velocity is, hence, given by 
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c(x) = 2 3 / 2 JQ(l - ^). (3.48) 



The averaged constraint (|2.3| ) governs, like for ferromagnets, the renormalized spin-wave 
gap Ai(x). For the correlation length at S 1 > 1 we find Eq. ( |3.11|) , with the exchange 
coupling renormalized according to Eq. ( |3.44j ): 



^ = ^[^(1--^)]. (3.49) 

As noted above, the SBMFA gives correctly only the exponent. The two-loop calculation 
of Chakravarty et alp as well as Monte Carlo simulations (Ref. [36] and references therein) 
and experimental datapH show that the prefactor C in this equation does not depend on 
the temperature up to small terms of the order of T/2iip s . At J' > J, we have 

e(T) = C7exp[?^^(l-^)]. (3.50) 

At sufficiently low temperatures the ratio £(x, T)/£(0, T) is small even at small dopant 
concentration x <C T/U <C 1. 

D. The phase transition in the quasi-2D antiferromagnets 

The decrease of the 2D correlation length with doping causes the rapid reduction of the 
Neel temperature in the quasi-2D antiferromagnets. Starting from the relation 1 

T N (x)^J ± e(T N ,x), (3.51) 

Eq. (|3l|D yields 

r»(*)=r„(o)[i- A(r= ' rw(0)) 1 , (3.52) 

where the Neel temperature for the undoped antiferromagnet is given by 

Ttv(0) w AtvJS 2 / \n(J/J ± ). (3.53) 
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Since x in Eq. (|3.52|) is multiplied by a large factor 1/Di(Tn(0)), the Neel temperature 

is supressed rapidly, when the dopant concentration increases. When J' > J, the Neel 

temperature decreases with the increase of doping as 

Ux 
T N (x)=T N (0)(l--—). (3.54) 

-LN{V) 

At a sufficiently large ratio J/J±, the Neel temperature extrapolates to zero at a doping 
concentrtion x which is much smaller than the percolation threshold. 

IV. FRUSTRATED PLAQUETTES 

In this section we consider the effect of a more complicated defect, i.e. a frustrated 
plaquette, on the properties of a 2D antiferro magnet. The plaquette is formed from 4 
frustrated bonds, which connect 4 neighboring spins. To begin with, let us summarize the 
results obtained in Appendix B in the harmonic spin-wave approximation. The T-matrix 
has a pole, when the frustrated coupling, J', reaches the critical value J c = 0.376 J. Like in 
the case of a single frustrated bond, the divergence of the T-matrix signals a local instability 
of the system. When J' > J c , the defect gives rise to a canted ground state. The threshold 
for instability of the perfectly aligned antiferromagnetic ground state is shifted to a weaker 
value, when the extra hole frustrates four bonds in a plaquette, rather than one bond. The 
reason is that in a plaquette two of the four bonds, connecting each spin with its neighbors, 
are frustrated. 

There is a remarkable difference between the effect of frustrated plaquettes on the spin- 
wave stiffness in the ferro- and antiferromagnets. In antiferromagnets the stiffness diverges 
when J' — > J~ . In ferromagnets the divergent term in the T-matrix scales as q A , and hence 
the stiffness passes smoothly the singular point J' = J c . Therefore we employ the SBMFA 
to consider frustrated plaquettes only in antiferromagnets. 

The interaction Hamiltonian we wish to treat is the following 

_ 4 

H^t = W p i Yl ( a i a m + aWm) + (Ai - A) ^ a\ a i- G 4 - 1 ) 

<lm> 1=1 
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Here the first sum < Im > runs over the frustrated bonds < 12 >, < 13>,<24>,<34>; 
Ai is the Lagrange multiplyer for the spins 1-4 at the corners of the plaquette. The energy 
W p i = J'Qpi + JQ, and Q p i is the correlation amplitude of the frustrated bonds. Evidently, 
Q p i is the same for all 4 frustrated bonds in the plaquette. 
The interaction energy matrix 



/ 



V 



\ 



(4.2) 



Ai-A W pl W pl 

W pl Ax -A W pl 
W pl A x -A W pl 

W pl Wpi Ai-A 

is diagonalized by the unitary transformation (B3). The Fourier-transformed interaction 
Hamiltonian (14.11) becomes 



V 



/ 



Him = J2 Wi J2 < 4f | yi{q) X yi{qi) \Af>. 

i=l qqi 



(4.3) 



Here the two-component column (row) vector | A^ > (< A^ |) is given by Eq. ( [3.13D 
w 1 =w 3 = JQijl p i - Wpi/2, w 2 = JQfipi - Wpi, w 4 = JQfj.pi = (Ai - A + 2W p i)/4, 



where 



yi(q) >-- 



( h(q) ^ 



rrii(q) 



<Vi(q)\=(im, rn*{q)) 



(4.4) 



f^q) = f 3 (q) = 1 - e - l ^ +q y\ mi {q) = -m 3 {q) = e~ iq * - e 



iq y 



(4.5a) 



h{q) = h(q) = 1 + e-^ +q y\ m 2 {q) = -m 4 {q) 



-iq* _ p-i<ly 



(4.5b) 



The equation for the one-defect T- matrix can be solved as before. The result is 

^/w - x A I Vi{q) ><yi{qi) I 



Ai(w) 



Here 



(4.6) 



Ai(w) = l-w&iiu), 



(4.7) 



2N 



where 



**(") = ^E < VU I £°(?» | y,($) > . (4.8) 



It is straightforward to calculate the functions ^i{ui) substituting expressions (|3.16|) , ( (4.4|) , 
( PI) into Eq. ( p|) . This yields 

* lM= * 3 H^E 1 ~T 1 T g " , (4.9a) 

M») = S -§Z S ™ 2 T S f"' , («b) 

, . . 8 JQ ^ 2 + cos 2 g x + cos 2 g„ + 4 cos g x cos q y . . 

^M = — E ^^ . (4.9c) 

Setting /i p ; = in Eq. (|4.6|) recovers in the quasiclassical limit the T- matrix expression 
derived in Appendix B, for the plaquette defect in the long range ordered antiferromagnet 
at zero temperature. 

Converting the sums over q into integrals over the first Brillouin zone in the reciprocal 
lattice, one finds numerically the values 

tfi(O) = -1.273/ JQ; tf 2 (0) = -0.727/ JQ. (4.10) 



Hence, 



A^O) = A 3 (0) = 0.124 + 0.876C P i + 1.273/v, (4.11a) 

A 2 (0) = Cpi + 0.727/ip,, (4.11b) 



where 



C P , = 1-0.727^ (4.12) 

Next we show that ( p i and /i p ;, like £ an d A 4 m the previous section, are small at low 
temperatures, if the frustration of the bonds in the defect plaquette is sufficiently strong, 
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i.e. if J' > J c , or | J' — J c | / J c <C 1. Hence, Ai(0) ^> A 2 (0), i.e. the first and third modes 
in Eq. ( |4.6| ) may be neglected when considering the effect of the frustrated plaquette. 

Given the expressions for the T-matrix and the functions \&i(co>), we can derive, as before, 
the equations for the functions ( p i(T) and /z p /(T). We find 

dun(uj){w 2 Im [ K \ . , n + wjm 1 K \ . . y v n } = 0, (4.13) 

3 A 2 (co>) A 4 (u;) 

Q P i = Q + ^ dun{u)Im 2V ;L 2 . v / ^A 4.14 

87T J-oo A 2 (c<j) 

Taking into account that at low frequencies, uj <C 4JS 1 , the functions \l/ 2 (u;) and ^4(0;) can 
be simplified as 

y 2 (u) = * 2 (0) + ^-g(u), y 4 (Lu) = 64JQg(uj), (4.15) 



we find from Eqs. ( |4.13|) and (|4.14| ): 



^V^^^ 1 ^ < 4 - 1& » 



t+ ^ + 1) = i/,m, (4.i6b) 



Q " (T) = Q j ' t(i+ ( l) + JC> ' (4 ' 1&) 



Here £=(</' — J c )/J c , j c = J c /J = 0.376, and 

1 r°° 
h(T) = — dxx 2 
n Jo 



1 r°° 

Tl JO 
XCOth^-M^,, -^-—r'\n~\' (I .y,.)^,- 1 }- (1.17) 



ZJLJX r| . /i p ; 1 + J c 2 1 l2 , -\2 4l— 1 

-77^1 Cp* + 7-7— x In— ■ ■■ ■ 

1 1 + J c 71 X 1 



Reasoning analogous to that given in the preceding section shows that the quantum 
corrections in Eqs. ( |4.16|) are of the order of (2S)~ 2 , when t > 0. In the large spin limit 
the functions (jl p i(T) and ( p i(T) are governed by equations, similar to Eqs. ( |3.33| ) for one 
frustrated bond in an antiferromagnet, 
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M& + y— -) - ^^(^7^1) (! + — y— ). ( 4 - 18a ) 



hUj^ + i) ri(r) 



t + l 16JQ 2 (1+ Jc )' 

When £ < and T/4JQ 2 <C t 2 <C 1 the solution of these equations is as follows: 



(4.18b) 



jc\t\ TUT) T 2 

Cpl = TTT c + iGJQ 2 (i+j c )\ty ^ ~ (i6J5 2 ) 2 1 1 1 <<c Cp '' (4,19) 



Tl{T) 
IQjc \t\JQ 2 



Q P i = Q(i- 77H ,;/ r0 2 )- ( 4 - 2 °) 



At positive £ and small temperatures T/4JQ 2 <C £ 2 , we find 

M 16JQ 2[ t(l+j c ) 2S* 3 (l + i)-S*Q 2 i J ' l ' ; 

^ ~ 16JS*2S**(l + t)-QH' { ] 

MQ) = <p, + ^= 't^S » (421c) 

l+j c t(l +j c ) 16JQ 2 

°* = FTT' 1 - le^'K + *> - 2^(1 + T- q^- 1} - (421d) 

In the zero temperature limit, Q p i = Q, if t < 0, and Q P 2 = Q/(t+ 1) < Q, iit > 0. Thus, at 
t = or J' = J c = 0.376J, the ground state changes from collinear to canted. It is argued in 
Appendix A that at large distances r from the defect plaquette, the correlation amplitude 
Qjtjt.g approaches to S as r~ 6 , i.e. the frustrated plaquette acts on the antiferromagnetic 
background at large distances as a quadrupole. 

We now average over the random distribution of the frustrated plaquettes, and obtain 
the renormalized spin-wave energy 

4 = wj[l - x(l + j c )Re—^-]. (4.22) 

31 



In the small q limit this becomes 

4=8J 2 Q 2 [l-*i±|]+A^ (4.23) 



where the gap A p i should be calculated self-consistently, as explained in Sec. prD} Eqs. ( |4.23| ) 



and ( |4.21c| ) yield the renormalized spin wave velocity for strongly frustrated plaquettes with 
t> 0, 

c(x) = 2 3 / 2 (l - ^f), (4.24) 

with 

8jQH(i+ 3c y 

Upl - l(T)(t+l) ■ (425) 

Note that the energy U p i is larger than the corresponding energy U (Eq. (|3.47|) ) for frustrated 
bonds. For strongly frustrated plaquettes, t > 0, at low temperatures T/4JQ 2 <C t 2 , the 
renormalized correlation length £ and Neel temperature T^(x) are thus given by 

£(7» = Cexp[^^(l - ^f)}, (4.26a) 

T N (x)=T N (0)[l-^-]. (4.26b) 



V. DISCUSSION. THEORY VS EXPERIMENT 

We have generalized the SBMFA for 2D doped magnets, which allowed us to study the 
effect of noninteracting arbitrary frustrated bonds on the spin-wave spectrum and on the 
2D magnetic correlation length at finite temperatures. We describe the defect bond by two 
local parameters, namely, the mean-field amplitude and the Lagrange multiplier. In a more 
rigorous treatment local mean-field parameters should be introduced also for the near-by 
bonds. However, even our simplified consideration is sufficient to describe correctly many 
peculiar properties of the doped 2D and quasi-2D antiferromagnets. 
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As noted in the Introduction, it is expected that in real cuprates the effective ferromag- 
netic coupling J' generated by the localized holes is larger than the coupling J ,113 i.e. J' 
exceeds the spin canting threshold J c (for frustrated bonds J c = J, while for plaquettes 
J c = 0.376J, in agreement with the results of the linear spin-wave theory, see Refs. |T6|.|r7| 
and Appendix B). In this case the renormalization of the antiferromagnet spin- wave velocity 
in the large spin limit is of the order of Ex/T, with E = U for frustrated bonds, and E = U v \ 
for frustrated plaquettes (see Eqs. (|3.48|) , fl3.47|) , Q4.24|) and (|4.25| )). Because of the large 



factor E/T the renormalization at low temperatures is large even at small doping level. 
It follows from Eqs. (|3.50|) and ( f4.26a|) that the effect of doping upon the correlation 



length, like upon the spin-wave velocity, increases with the decrease of the temperature. 
The ratio £(x, T)/£(0,T) is of the order of unity if 2np s E/T 2 <C x~ l , and is exponentially 
small when {2np s Ex) 1 ^ 2 ^> T ^> Ex. The rapid decrease of the 2D correlation length 
strongly reduces the quasi-2D Neel temperature (see Eqs. ( |3.54| ) and (f4.26b| ). The above 



results hold also for frustrated bonds in ferromagnets (Eqs. ( |2.58| ), ( |2.59| ), ( [2.49 b| )). In 



contrast, frustrated plaquettes have only weak effects on the properties of the ferromagnet, 
the renormalization of the spin wave stiffness being of the order of x and temperature 
independent (Appendix B). 

The effect of frustrated bonds is large also at J' < J c , provided | t | = (J c — J')/ J c <C 1. 
However, the temperature dependence is weaker than at J' > J c . The renormalization of 
the spin wave spectrum in the large spin limit is proportional to x/T 1 / 2 if T 3> AJS 2 t 2 , and 
it tends to a constant value x/ \ t |> x when T -> (Eqs. (pT57|) , (|2T46|) , (|2~47a|) , (|2~48al) ). 



The strong effect of frustrating defects on the spin-wave stiffness and on the correlation 
length is related to the local instability leading to the spin canting. Indeed, the spin- 
wave scattering amplitude at T = uj = diverges at the local stability threshold J' = J c . 
At J' > J c the defects act as dipoles (bonds) or quadrupoles (plaquettes). Hence, the 
disturbance of the background decays as a power of the distance from the defect rather than 
exponentially. This means that the scattering amplitude has a zero frequency pole at any 
J' > J C .E3 At finite temperatures the scattering amplitude is finite but large and increases 
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with the decrease of the temperature. We believe that these arguments hold for any type 
of frustrating defects, and therefore the above results are valid qualitatively not only for 
frustrated bonds and plaquettes. 

We also calculated the quantum corrections to the local mean-field amplitude in anti- 
ferromagnets in the large spin limit. The corrections are relatively large, ~ (S 2 \nS)~ 1 / 3 , 
near the threshold (Eq. ( |3.37|) ), however at stronger frustration (J' > J c ) the corrections are 
unusually small, ~ (2S') _2 ln~ 1 2S (Eqs. (|3.38| )) and (|3.39|) ). This result gives evidence in 
favour of the classical model for strongly frustrated bonds in 2D antiferromagnets, proposed 



in Ref. 12 



It follows from our theory that weakly frustrating and nonfrustrating impurities act quite 
differently. The renormalization of the spin stiffness is of the order of x <C 1, and does not 
depend on the temperature. Thus, the correlation length is given by the same Eq. (1.1), as 
for undoped samples, with a concentration dependent constant A. Similar results have been 
obtained in Ref. |37] by means of the quantum Monte-Carlo method. 

Let us compare our results with the experimental findings&aMII^SE] for the temperature 
and doping concentration dependences of the correlation length in quasi-2D antiferromag- 
nets. It has been shown in these papers that one should distinguish between two kinds of 
impurities. Dopants which introduce holes and, presumably, frustrate the Cu-Cu bonds, 
supress strongly the correlation length at low temperatures and reduce rapidly the Neel 
temperature. In contrast, dopants which introduce vacancies or excess electrons do not 
change the temperature dependence of the correlation length and simply reduce the stiffness 
in accord with percolation. The reduction of T/v is much more gradual than that in hole 
doped samples. The results, presented here, account qualitatively for all these properties of 
doped antiferromagnets.E3 

A remarkable difference between the magnetic properties of the electron- and hole-doped 



materials was observed in Ref. ^9[ In electron doped Nd2- x Ce x CuO^ the reduction of T/v 
and p s is of the same order of magnitude. However, in the hole doped LaiCuO^+x the 
correlation length (and, hence, p s ) at high temperatures (T « 500-K" ^> T/v(x)) is almost 
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unchanged by doping, while Tn essentially decreases. Our results permit to explain this 
puzzle. In electron doped samples the effect of doping does not depend on the temperature, 
and the decrease of p s , as well as of £ and T/v, is of the order of x. In hole doped samples 
the change of T/v is larger than that of £(T ^> T/v), since the effect of doping grows with the 
decrease of the temperature. 

We now show that the experimental results confirm the temperature dependence of the 
renormalized correlation length, derived in the preceding sections. It follows from Eqs. (|3.50|) 



and (}4.26a|) that in antiferromagnets, doped with strongly frustrating impurities, the quan- 
tity Tm(£/C) should be a linear function of x/T: 

T\J-^=2n Ps (l-^), (5.1) 



where E, as before, stands for U or U v \. In Refs. p|,|9|JT0| the correlation length was measured 
for three doped samples LaiCuO±+ x , with T^[x) = 245 K, 190 K and 90 K. The precise hole 
concentrations for these samples are not given in the references cited above. It is known, 
however ,@'e21 that the hole concentration varies approximately linearly with 325 — Tjv(x). 
Therefore one has an estimate for their ratios, x(90) : x(190) : x(245) = 1 : 0.59 : 0.34. 
The quantity C can be determined from the measurements of £(T, 0) in an undoped sample 2 
( wee neglect its week temperature dependence, see Sec. | [IQ) , and the stiffness constant 



is 2np s ~ 150 meV.H Using these data, we plotted in Fig. 1 the experimental values of 
Tln(£/C) vs x/T. It is seen that the experimental points for all three samples fall on a 
single straight line, as expected from Eq. (|5.1|) . 

According to Ref. |2|, the linear function T/v(x) extrapolates to zero at x ~ 2%. Thus, 
we can estimate the hole concentration x(90) ~ 1.4%. Then, from the slope of the straight 
line in Fig. 1 we estimate the energy E ~ 500 meV. The theoretical values, which follow at 
t > 1 from Eqs. ( gg7p , ( ggg ) and the relation! 2tt Ps = 1.15 J, are: U = 110 meV, U pl = 
410 meV. The last number is close to the above experimental value. Note that the extra 
hole in lanthanum cuprate is localized on a region of the order of two lattice constants,^ 
and hence the real defect is more complicated than the above models. What is more, our 
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study is based on the assumption that the impurity holes are localized, which is, apparently, 
violated at high temperatures.EJc3 Nevertheless, the above analysis shows that the theory 
of strongly frustrating defects explains the experimental findings in samples with Neel order 
at least semiquantitatively. 

Keimer et alB also measured the temperature dependence of the correlation length £(T, x) 
of Sr doped samples (x= 0.02; 0.03; 0.04), without Neel order. Eq. ( |5.1|) fails to account for 
these experimental data even at high temperatures, when Ex/T < 1. Keimer et alB showed 
that in this case the correlation length obeys the empirical relation 

r\x,T) = r\x,o) + r l (o,T), (5.2) 

where £(0, T) is the correlation length of the carrier-free sample, and £(x, 0) is the measured 
correlation length at low temperatures (£(x, 0) = 150, 65 and 45 A for x = 0.02, 0.03 
and 0.04 respectively). They also attemped to fit Eq. ( |5.2|) to the data for the above 
less doped samples, which exhibit 3D long-range order, £(x, 0) being a fitting parameter 
(£(x, 0) = 140 and 275 A for the T^ = 90 and 190 K samples respectively). The fit is good 
only at relatively high temperatures. Figure 1 demonstrates that Eq. ( p.l| ) describes the 
experimental data better than Eq. Q5.2|), especially for the T N =90 K sample, although one 
needs only one parameter (E) to fit the data for all samples. This implies, first, that £(x, 0) 
for the less doped samples is larger than the measured £(x, T), namely, £(x, 0) > 400 A for 
the T/v = 90 K sample. One cannot rule out the possibility that £(x, 0) is infinite for the 
concentrations x smaller than a certain (small) critical concentration x c . This standpoint is 
somewhat supported by the numerical simulations in the zero-temperature limit .c3a It has 
been shown that the frustrated bonds destroy the 2D long-range order in the Heisenberg 
antiferromagnet (£(x, 0) is, perhaps, finite at any x"), while in the XY model the correlations 
do not decay exponentionallay, but rather as a power of the distance.^ One can, thus, suggest 
that the small XY-anisotropy, which exists in the real cuprates, stabilizes the 2D, T = 
long-range order in doped samples at some small, but finite doping concentration. 

The finite value of £(x, 0) in the samples with high doping level also affects the renor- 
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malization of the spin-stiffness. We showed in Sections [ HIC| and [TV] that the stiffness renor- 



malization given by Eqs. ( |3.46| ) and ( |4.23| ) holds only in the small q limit, when 

q2 <<: ehLs/^y (5 - 3) 

while at larger wave vectors the spin-stiffness is the same as in the undoped samples. On 
the other hand, the spin waves exist only within a region of size £(x, T), i.e. 

g»r 1 (^,T) (5.4) 



(See Ref. ^6] and references therein). In slightly doped samples, the inequalities (|5.3| ) and 
( |5.4| ) are always fulfilled, since the inverse correlation length is exponentially small at low 
temperatures. In strongly doped samples, a new temperature independent length scale 
£(x, 0) appears. Thus, the region of well-defined spin waves is restricted to distances smaller 
than £(x, 0), while according to Eq. (|5.3j ) only large wavelength spin waves are renormalized. 
At sufficiently high doping level, i.e. small £(x, 0), the inequalities ( |5.3| ) and ( |5.4| ) are 
violated, and hence the spin-stiffness is not renormalized. This reasoning explains why the 
correlation length for the x = 0.02 sample is larger than for the T/v = 90 K sample (compare 
Figs. 8b and 10 in Ref. 0). 

Substituting an estimated value of the energy E p» 500 meV into Eq. ( |3.54|) for the 



Neel temperature, we find T N (x) = T N (0)(1 — bx), where b « 20 per 1%. The coefficient 
b is smaller than the experimental value by a factor 2.5 — 3. There are several reasons 
for this discrepancy. First, Eq. Q3.54J ) holds only at small x when bx <C 1. Second, all 
our numerical estimates are based on the above approximate linear relation between the 
function 325 — T^{x) and x. Finally, the mean- field equation ( |3.51| ) which relates the Neel 
temperature to the 2D correlation length should be improved.^ 

As noted above, the holes in La2CuO^ +x are delocalized at temperatures T > Tjy. The 
overall agreement of our theory, based on the frustration model, with the experimental 
results for hole doped lanthanum cuprate supports the standpoint that even delocalized holes 
frustrate the antiferromagnet. It would be difficult to understand the striking difference 
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between hole and electron (vacancy) doping, as well as the temperature dependence of 
the spin stiffness in hole doped samples, without this assumption. This conclusion agrees 
with the results of our analysis^ of the reentrant temperature dependence of the sublattice 
magnetization in oxygen doped cuprates. We gave evidence there that delocalized holes 
frustrate the antifferomagnet, though less than localized holes. New efforts are needed in 
order to elucidate this question. 



Finally, we would like to note that, as has been shown in Sections |TJ and [TV], the spin- 
wave velocity c should also change drastically with hole doping. Even though we considered 
only 2D antiferromagnets, we expect that the same renormalization takes place in quasi- 
2D antiferromagnets at T < T/v for spin waves with wave vectors in the range J±/J <C 
q 2 <C T/E\n(4:JS/ujg). Here the first inequality guarantees that the interplane coupling 
does not affect the spectrum, while the second one coincides with Eq. ( |5.3|) . Thus, we may 
conjecture that the renormalization of c in hole doped quasi-2D antiferromagnets is of the 
order of JS 2 /T in the large spin limit. The same conclusion follows from a direct calculation 
of the spin wave scattering on free dipoles at T <^C T/v .^3 Unfortunately, the experimental 
information on this subject is very poor. We are aware only of the Ref. fi5| , fi6| , where it 
is shown that a small amount of holes in the CuOi planes of YBa2Cu^O^ +x renormalizes 
strongly the spin-wave velocity. Neither the temperature nor the concentration dependences 
of the spin-wave velocity has been studied. We hope that our results will stimulate new 
experiments in this field. 
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APPENDIX A: THE CORRELATION FUNCTION Q., 



f,f+5 



In this appendix we calculate the correlation of the spins at the ends of a bond, 5, at a 
large distance from the frustrated bond, a. 

1. Ferromagnet 

We start from the expression 

Qr,?+S = < a W+s > 

= — / dun(u) — ^ImG&fau^*-*****, (Al) 

7T J— oo iV .— 

and subtract the constraint equations for the spins at the ends of the bond S, 

/+oo 
d u n{u) £ e^-^Im[G(q, q x] u) - G°(q, w)%] 

/+oo 
dun{u)Y^ qi '^ r+s) Im[G{q,q x] u) - G°(g»%] = 0. (A2) 

-oo -« 

We find 

1 r+oo I 

<W = 5 - — J dun^—ImYe^-^lG^q^uj) - G°(q,u;)8 q ^} 

x(l-e-^)(e^-l). (A3) 

Eqs. (K3), (gTg), and (|Q5Q yield 



1 /-+00 I _ 

Q g =S + - dun(u)Im-—-[Y / G°(q,u;)(q-d)(q-S)e i ^ (A4) 

Z7T J — oo ivx^CJJ 

where we have taken into account that at large r only small wave- vectors are important, and 
expand the exponents in powers of (q- a), (q- 5). The term proportional to v-i = \Xf gives at 
low temperatures a negligible small contribution to Q ? ? + $. At small frequencies 10 <C JS/r 2 , 
the integration over q in Eq. (|A"4]) is easily performed, and we find 



Qr,f + s =3- J^f F ^Z f {T). (A5) 
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Here 

1 r- * 2 



F ^ -- zh&'S-iktf' %(f- 5 )] 2 > ( A6 ) 



ifr £ <y"-i tyt A 



1 f+°° 1 

Z f (T) = — doon^Im—— (A7) 

ZTlJ-oo DAuj) 



It follows from Eqs. ( |2.48a| ), ( |2.49b| ), and ( |A7|) that in the zero temperature limit 



t + 1 

Thus, we finally obtain that Qf^ + g = Q, when t < 0, and 



Z f (T) = for t < 0, 

Z^T) = =1^-1 for t > 0. (A8) 



^r 5 ! 1 -^! (A9) 



when t > 0. The r dependence of Q ?if+ ^ given by Eqs. ( |X9| ) and ( |A6|) holds in the region 
1 < r < (AJS 2 /Ty/ 2 . 

2. Antiferromagnet 

We first consider a frustrated bond. In this case the function Qf^ can be written as 

J r+oo ] ^ ^ 

^r+5 = Q / dujn(uj)Im— J2 < Xi I [ G (<?> ?i! w) - G°(?»%] | Xt > 

7T J—oo iV —. 

991 

x e <(7-ft)*+^. (A10) 

Taking into account the constraints for the spins at the ends of the frustrated bonds, we 
rewrite this equation as 



W f +OD 1 



1 ™ «.Jf ,An,. , 

X 



< xt I G°(g» | Vl (q) >< ^(gl) | G°(gi,u;) | Xt > (1 + e «t-90-*)] e *t**K (All) 
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We expand the expression within the square brackets in Eq. ([All ) in powers of (q ■ a) and 
(q ■ 5), and transform it into a sum of three terms 

-[±J2e iq "G° n (M(q- a)(q- <T)] 2 + {| ^ e^[G° n (q, ") " G? 2 (g»]} 2 

q q 

M 
+^E[^i(g^)-G? 2 (g,c,)]e^5:te-«)G? 1 (gi^)e^-. (A12) 

q qi 

The first of these terms, when integrated over q, transforms in the zero frequency limit into 
the above function (|A6|) . The last two terms have no analogy in the case of ferromagnets. 
Their contribution to Q?jt + g can be important at small r. However, at large distances, one 
has 



jdq\G\M - chW**" j d *\jr ei * r ~= °- ( A13 ) 



,1 -JT9 
q 

Thus the behavior of the function Q^p,^ at large distances in antiferromagnets is the same 

as in ferromagnets: S — Q?? + $ oc r~ 4 in the region 1 <C r < (4J5' 2 /T) 1//2 . Analogous 

reasoning holds also in the case of frustrated plaquettes. Then, instead of the first term in 

Eq. (|ATJ) we find 

[Y,e^GUm-a){q-b){q-5)]\ 
q 

where a and b are unit vectors along the plaquette sides. Hence, at large distances the 
function S — Q^f + g decays as r~ 6 . 

APPENDIX B: THE FRUSTRATED PLAQUETTE AT ZERO TEMPERATURE. 
HOLSTEIN-PRIMAKOFF REPRESENTATION 

In this Appendix we show that the effect of the frustrated plaquettes on the spin-wave 
spectrum is strikingly different in ferro- and antiferromagnets. While in antiferromagnets 
the renormalization diverges when the frustrated coupling J' tends to the local stability 
threshould, J c , in ferromagnets the renormalization is finite in this limit, and, hence, is 
qualitatively the same at J' < J c and J' > J c . 
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1. Ferromagnet 

We suppose that the strength of the frustrated bond J' is smaller than the threshold 
value J c for local instability (the value of J c is obtained below). Then all the spins are 
parallel, and the Hamiltonian, which describes the interaction of the spin waves with the 
spins at the corners of the frustrated plaquette, is 



Him = (J + J')S J2 K ~ a ln)( a m ~ &l 



(Bl) 



<lm> 



where < Im > is summed over the bonds in the plaquette. It follows from ( |B1| ) that the 
interaction matrix is 

'-2 1 1 



V = (J+J')S 



V 



1-201 
10-21 
11-2 



(B2) 



It is easy to verify that the matrix P = UVU 1 , where 



'i-i i -^ 



U 



1-1-1 1 

11-1-1 

1111 



(B3) 



is diagonal, and its elements are Pi = P 3 = —(J + J')S/2,P 2 = —(J + J')S,Pi = 0. The 
Fourier transform of the Hamiltonian ( [BID is given by 



Hint = Yl P i^ Yl X i(0) X i(qiWqa 



1=1 qq 1 



q^qii 



(B4) 



with 



Xl {q) = (1 - e - i «' s )(l + e~^), 



(B5a) 



x 2 (q) = (1 - e-^)(l - e- iqb ) 



(B5b) 
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x 3 (q) = (l + e-^ a )(l - e- l * b ). (B5c) 

The solution of the T-matrix equation can be obtained in the usual way, 

x<(g)xj(gi 



where 



T(g,g 1 ^) = gP li _ p ^ M , (B6) 



q 



It follows from Eqs. ( |B7| ) and ( |B5| ) that 

= , x * / \ , f d 2 q (1 + cos Oa;)(l — cos g„) ,„„ , 

$! w = $3 w = 4 / — -^ ^ ^, (B8a) 

<& 2 (0)=4/ , /? 9 (1 " COSfe)(1 " COS ^ ) . (B8b) 

Thus, the function T(g, g; a;) is given in the limit of smal q and uj by 

7Y--rrt 4Pl 1 2 P2( & q2 v raa\ 

T(g ' q] 0) = " i-px^o) " i-p 2 * 2 (or (B9) 

The second term in the r.h.s. of this equation is of higher order in small q, and can be 
neglected at small P 2 . However, it appears that with the increase of P 2 it is the denominator 
of this term which approaches to zero first, and, hence, determines the threshold for the local 
instability 

The integration over the first Brillouin zone of the reciprocal lattice yields 

$ 2 (0) = -0.727(J + J')/JS, $x(0) = -1.363/ JS. 

Hence Pi$i(0) = 0.682 (J + J')/ J < P 2 $ 2 (0) = 0.727 (J + J')/ J, and the threshold value 
of J' should be found from the equation 

1-P 2 $ 2 (0) = 0. 

Thus J c = 0.376J, which exactly coincides with the value of J c , found in Sec. [TV] in the 
Schwinger boson representation for the plaquette defect in an antiferro magnet. 
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Note that the renormalized spin-wave spectrum 

e g - =uj q + xT(q, q) = u ? [l - -^^ ,., x\ 



(BIO) 



J(J C - 0.806J') 

is nonsingular at J' = J c . Therefore the renormalization of the spin-stiffness in the ferro- 
magnet with frustrated plaquettes, unlike ferromagnets with frustrated bonds, is small, of 



the order of x, at both J' < J c and J' > J c . 



2. Antiferromagnet 

The Hamiltonian, which describes the interaction of the spin waves with the spins at the 
corners of the frustrated plaquette, when J' < J c , is 



Hi nt = ~(J + J')S Yl ( a l b ™ + <4&™ + a\di + b^bm), 
<lm> 



(Bll) 



where ai(b m ) are the Holstein-Primakoff operators for the spins in the sublattice A (B). The 
interaction matrix, like in ferromagnets, can be diagonalized by the transformation 
Then the Fourier transform of the Hamiltonian is given by 



Here 



3 1 

H in t = E P i^ £) < C 9" I **($ >< Y i&) I C 9"l > • 
t=l g§i 



Cq*> = 






< Cn 



(4> &-*) 



(B12) 



(B13) 



where 



*(g) >= 






< y,(g) |= (a*(q) (3m) 



(B14) 



«i(g) =a 3 (q) = l-e 



-»(9a:+9 W ) 



Pi® = -P3(q)=e- Ub -e**>, 



(B15a) 



0:2(9) = 1 + e 



*(#■:+?!/) 



p2{q) = e~ iq '+e 



«/;, 



(B15b) 
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The next calculations are straightforward, provided we introduce the Green's function matrix 
by Eq. (|3.15| ) with the vectors < cg\ (\ cg>) substituted for < A^ | (| A^ >). 
The unperturbed Green's functions are 



G°(q,cu) 



u 2 - U)\ 



(B16) 



US+LJ -US-fq 

y -4JS 7g - 4 JS - UJ j 

Like in the ferromagnet, only the mode i = 2 is relevant at J' near the threshold J c = 0.376 J. 
Therefore the single-defect T-matrix can be written as 

T( q , qi ; u>) = r^RW) ' (B17) 

where 

_ 8 JS v sin 2 q x + sin 2 q y 
R^)-^f\ u2 _ 2 • (B18) 

The function R(u) coincides up to the quantum corrections with the function ^ 2 { UJ ) from 
Eq. flOEp . Thus i?(0) = 0.727(1 + J /J'). The Eqs. (|BT7|) , ( fUg) , and (pT5|) together with 
the relation,c3 

x 
e g - = w ? + — [4JS , (Tn(g, <f) + T 22 (q, q)) + (jJ^T u (q, q) - T 22 (q, q)) 

ZUq- 

- AJS lg iT 12 (q, q) + T 21 (q, q))} (B19) 

yield the renormalized spin-wave spectrum, 

The remarkable difference between the effect of the frustrated plaquettes on the spin-wave 
spectrum in the ferro- and antiferromagnets is now evident. 
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FIGURES 

FIG. 1. (T/2TTp s ) ln(£/C) as a function of 2iTp s x/x(90)T for three samples 
LaiCuO^+x with T/v= 245 K (x), 190 K (o), and 90 K (•). The points are obtained from 
the experimental data (Refs. |2||],|II]), as discussed in the text. The straight solid line shows 
the fit to Eq. (|5.1|). The dotted and dashed lines represent fits to Eq. ( p.2| ) for the T N = 90 
and 190 K samples (£(z, 0) = 140 and 275 A, Ref. g). 
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